For requsting MC pp#2350
Conversation
|
REQUEST FOR PRODUCTION RELEASES: This will add The following labels are available |
jackal1-66
left a comment
jackal1-66
left a comment
There was a problem hiding this comment.
if (child0 != nullptr && child1 != nullptr) {
// check for parent-child relations
auto pdg0 = child0->GetPdgCode();
auto pdg1 = child1->GetPdgCode();
std::cout << "First and last children of parent " << checkPdgSignal << " are PDG0: " << pdg0 << " PDG1: " << pdg1 << "\n";
if (std::abs(pdg0) == checkPdgDecay && std::abs(pdg1) == checkPdgDecay && pdg0 == -pdg1) {
nLeptonPairs++;
if (child0->getToBeDone() && child1->getToBeDone()) {
nLeptonPairsToBeDone++;
}
}
Your test gets stuck here which could mean either that your generator is broken and does not generate the correct pairs, or something is fishy in the test itself
| << "#signal (prompt Jpsi): " << nSignalJpsi << "; within acceptance " << rapiditymin << " < y < " << rapiditymax << " : " << nSignalJpsiWithinAcc << "\n" | ||
| << "#signal (prompt Psi(2S)): " << nSignalPsi2S << "; within acceptance " << rapiditymin << " < y < " << rapiditymax << " : " << nSignalPsi2SWithinAcc << "\n" | ||
| << "#lepton pairs: " << nLeptonPairs << "\n" | ||
| << "#lepton pairs to be done: " << nLeptonPairs << "\n"; |
There was a problem hiding this comment.
This should be nLeptonPairsToBeDone
| << "#antileptons: " << nAntileptons << "\n" | ||
| << "#signal: " << nSignal << "\n" | ||
| << "#lepton pairs: " << nLeptonPairs << "\n" | ||
| << "#lepton pairs to be done: " << nLeptonPairs << "\n"; |
There was a problem hiding this comment.
This should be nLeptonPairsToBeDone
| auto pdg0 = child0->GetPdgCode(); | ||
| auto pdg1 = child1->GetPdgCode(); | ||
| std::cout << "First and last children of parent " << checkPdgSignal << " are PDG0: " << pdg0 << " PDG1: " << pdg1 << "\n"; | ||
| if (std::abs(pdg0) == checkPdgDecay && std::abs(pdg1) == checkPdgDecay && pdg0 == -pdg1) { |
There was a problem hiding this comment.
Equivalent to pdg0 == -pdg1 && std::abs(pdg0) == checkPdgDecay
| auto pdg0 = child0->GetPdgCode(); | ||
| auto pdg1 = child1->GetPdgCode(); | ||
| std::cout << "First and last children of parent " << checkPdgSignal << " are PDG0: " << pdg0 << " PDG1: " << pdg1 << "\n"; | ||
| if (std::abs(pdg0) == checkPdgDecay && std::abs(pdg1) == checkPdgDecay && pdg0 == -pdg1) |
There was a problem hiding this comment.
Equivalent to pdg0 == -pdg1 && std::abs(pdg0) == checkPdgDecay
2 participants
No description provided.